function w_mult = f_w(w)
    global alpha
    w_mult = (w  ).^(-alpha);
end